Rome III Criteria Capture Higher Irritable Bowel Syndrome SNP-Heritability and Highlight a Novel Genetic Link with Cardiovascular Traits

Background & Aims Irritable bowel syndrome (IBS) shows genetic predisposition, and large-scale genome-wide association studies (GWAS) are emerging, based on heterogeneous disease definitions. We investigated the genetic architecture of IBS defined according to gold standard Rome Criteria. Methods We conducted GWAS meta-analyses of Rome III IBS and its subtypes in 24,735 IBS cases and 77,149 asymptomatic control subjects from 2 independent European cohorts (UK Biobank and Lifelines). Single-nucleotide polymorphism (SNP)-based heritability (h2SNP) and genetic correlations (rg) with other traits were calculated. IBS risk loci were functionally annotated to identify candidate genes. Sensitivity and conditional analyses were conducted to assess impact of confounders. Polygenic risk scores were computed and tested in independent datasets. Results Rome III IBS showed significant SNP-heritability (up to 13%) and similar genetic architecture across subtypes, including those with manifestations at the opposite ends of the symptom spectrum (rg = 0.48 between IBS-D and IBS-C). Genetic correlations with other traits highlighted commonalities with family history of heart disease and hypertension, coronary artery disease, and angina pectoris (rg = 0.20–0.45), among others. Four independent GWAS signals (P < 5×10-8) were detected, including 2 novel loci for IBS (rs2035380) and IBS-mixed (rs2048419) that had been previously associated with hypertension and coronary artery disease. Functional annotation of GWAS risk loci revealed genes implicated in circadian rhythm (BMAL1), intestinal barrier (CLDN23), immunomodulation (MFHAS1), and the cyclic adenosine monophosphate pathway (ADCY2). Polygenic risk scores allowed the identification of individuals at increased risk of IBS (odds ratio, 1.34; P = 1.1×10-3). Conclusions Rome III Criteria capture higher SNP-heritability than previously estimated for IBS. The identified link between IBS and cardiovascular traits may contribute to the delineation of alternative therapeutic strategies, warranting further investigation.

BACKGROUND & AIMS: Irritable bowel syndrome (IBS) shows genetic predisposition, and large-scale genome-wide association studies (GWAS) are emerging, based on heterogeneous disease definitions.We investigated the genetic architecture of IBS defined according to gold standard Rome Criteria.

METHODS:
We conducted GWAS meta-analyses of Rome III IBS and its subtypes in 24,735 IBS cases and 77,149 asymptomatic control subjects from 2 independent European cohorts (UK Biobank and Lifelines).Single-nucleotide polymorphism (SNP)based heritability (h 2 SNP ) and genetic correlations (r g ) with other traits were calculated.IBS risk loci were functionally annotated to identify candidate genes.Sensitivity and conditional analyses were conducted to assess impact of confounders.Polygenic risk scores were computed and tested in independent datasets.RESULTS: Rome III IBS showed significant SNP-heritability (up to 13%) and similar genetic architecture across subtypes, including those with manifestations at the opposite ends of the symptom spectrum (r g ¼ 0.48 between IBS-D and IBS-C).Genetic correlations with other traits highlighted commonalities with family history of heart disease and hypertension, coronary artery disease, and angina pectoris (r g ¼ 0.20-0.45),among others.Four independent GWAS signals (P < 5Â10 -8 ) were detected, including 2 novel loci for IBS (rs2035380) and IBSmixed (rs2048419) that had been previously associated with hypertension and coronary artery disease.Functional annotation of GWAS risk loci revealed genes implicated in circadian rhythm (BMAL1), intestinal barrier (CLDN23), immunomodulation (MFHAS1), and the cyclic adenosine monophosphate pathway (ADCY2).Polygenic risk scores allowed the identification of individuals at increased risk of IBS (odds ratio, 1.34; P ¼ 1.1Â10 -3 ).
I rritable bowel syndrome (IBS) is a common disorder of gut-brain interaction. 1 It affects women more often than men and has a prevalence of 3%-11% in European populations, based on the consensus Rome criteria from the Rome Foundation. 2 These criteria define IBS as a recurrent abdominal pain or discomfort associated with abnormal stool frequency and consistency, and classify IBS subtypes as constipation (IBS-C), diarrhea (IBS-D), or mixed (IBS-M) based on the predominant stool type. 3The latest Rome IV criteria are the most stringent, focusing only on abdominal pain (rather than also discomfort, as in Rome III) and an increased frequency of symptoms (1 day/week vs 3 days/ month in Rome III), thus possibly capturing a more severe phenotype compared with earlier versions. 4,5xisting treatments for IBS are mostly directed toward the amelioration of symptoms, but are not effective in all patients. 6The contribution of dietary and psychological factors, gut dysbiosis, prior gastrointestinal (GI) infection, and genetic predisposition is recognized. 1 A heritable component of IBS has been demonstrated in twin and family studies, although with varying degrees of estimates depending on the adopted IBS definition. 7The highest heritability was estimated in the Australian Twin study, with 57.9% (95% confidence interval, 40.6-75.9) of the variance in the reporting of IBS symptoms being attributed to additive genetic effects. 8Whether these genetic factors differently predispose to various IBS forms and/or subtypes is unclear, which hampers personalized approaches to treat patients with distinct clinical presentations. 6nterest in the genetic architecture of IBS has been growing, because it holds potential to reveal actionable pathways and biomarkers. 9This can be sought, for instance, by analogy with other diseases showing similar pathophysiologic mechanisms, including those accounted for by shared genetic factors.5][16][17] However, these studies used heterogeneous IBS definitions that show poor overlap and inconsistent findings across definitions.For instance, the largest GWAS meta-analysis included results from 53,400 IBS cases, collectively identified through hospital-inpatient records, self-reported diagnoses, and Rome criteria: among those included from the UK Biobank (UKBB), less than 1% (n ¼ 340/40,548) satisfied all IBS definitions simultaneously. 15The applicability and clinical relevance of genetic findings obtained using pooled IBS definitions are yet to be determined, nor have different IBS manifestations been adequately explored based on stratification of patients in different Rome subtype groups.
Here, we investigate the genetic underpinnings of IBS and its subtypes, based on available genotype and Rome III questionnaire data from 2 large European population-based biobanks (UKBB and Lifelines [LL]), totaling 24,735 cases and 77,149 asymptomatic control subjects.

Characteristics of Rome III IBS Patients
Rome III IBS patients and asymptomatic control subjects were defined based on digestive questionnaire data available for UKBB and LL participants (see Methods, Figure 1).In total, we identified 24,735 Rome III IBS patients (22,745  UKBB patients and 1990 LL patients) and 77,149 asymptomatic control subjects (66,631 from UKBB and 10,518 from LL).The demographics of study participants are summarized in Table 1.In line with previous studies, IBS-M was the most common subtype both in UKBB (54% of patients with IBS) and LL (40%).Patients with IBS were predominantly women, with the IBS-C group showing the most pronounced female predominance (84.0% of IBS-C cases were female in both UKBB and LL). 18,19ealth-related data available in UKBB allowed other diseases and conditions to be studied in IBS cases and control subjects, based on hospital-inpatient records in relation to diseases (International Classification of Diseases, 10th version [ICD10], Chapters I-XIV) and laboratory findings (ICD10 Chapter XVIII), across 22,344,000 data-points (250 data-fields covering 89,376 individuals) (Methods).Age-and sex-adjusted regression analysis revealed patients with IBS (all subtypes) to suffer from a significantly higher number of comorbidities compared with control subjects (P 5.1Â10 -161 ) (Table 2).At both the chapter (domain) and 3-digit (trait-specific) ICD10 levels, strongest associations were detected for GI (chapter XI diseases of the digestive system) and psychiatric (chapter V mental and behavioral disorders) disorders, whose risk was more than doubled in Rome III IBS versus control subjects (odds ratio [OR], 2.0-2.7;P 7.7Â10 -39 ) (Table 3), particularly for F41 other  anxiety disorders (OR, 3.3-3.9;P 4.5Â10 -26 ) and K30 dyspepsia (OR, 2.9-3.8;P 1.0Â10 -44 ) (Figure 2, Supplementary Table 1).As a novel observation, traits from the chapter IX diseases of the circulatory system showed increased prevalence among patients with IBS (Figure 2, Supplementary Table 1).This included I10 essential (primary) hypertension, which was associated with all IBS subtypes (OR, 1.4-1.6;P 1.3Â10 -12 ), more pronouncedly IBS-M (Figure 2).Other differences, also related to diseases of the circulatory system, were observed across subtypes including I20 angina pectoris and I25 chronic ischemic heart disease, which were more common in IBS-C patients (OR, 2.2 [P ¼ 1.2Â10 -15 ] and OR, 1.8 [P ¼ 2.1Â10 -11 ], respectively) (Figure 2, Supplementary Table 1).

Rome III IBS GWAS Meta-Analyses: Heritability and Genetic Correlation with Other Traits
Individual Rome III IBS GWAS analyses were carried out in UKBB and LL, adopting a common analytic pipeline for quality control (per-sample and per-marker), imputation, and association testing adjusting for age, sex, and principal components (PCs) (Methods).UKBB and LL GWAS summary statistics for 6,311,313 high-quality single-nucleotide polymorphism (SNP) markers were included in GWAS metaanalyses spanning, respectively, 24,735 IBS, 13,132 IBS-M, 6481 IBS-D, and 4330 IBS-C cases and a common set of 77,149 asymptomatic control subjects (Table 1).GWAS meta-analysis results showed no population stratification (l¼1.03-1.09,linkage disequilibrium score regression [LDSC], intercept¼0.97-1.00).Based on LDSC analyses (Methods), all Rome III IBS subtypes showed detectable and significant SNP-based heritability (h 2 SNP ), estimated at 10.3% (P ¼ 7.1Â10 -25 ) for IBS, 13.4% (P ¼ 3.2Â10 -15 ) for IBS-M, 12.4% (P ¼ 1.7Â10 -8 ) for IBS-D, and 11.7% (P ¼ 1.6Â10 -4 ) for IBS-C.Considerable genetic overlap was observed across all traits, with almost complete correlation between IBS and IBS-M (r g ¼0.98; P < 1.0Â10 -300 ) (Figure 3, Supplementary Table 2).Of note, although weaker than with other subtypes, positive genetic correlation was also observed between traits at the opposite ends of the stool consistency spectrum, such as IBS-C and IBS-D (r g ¼0.48; The genetic architecture of IBS was then compared with that of other traits and diseases via LDSC analyses implemented in the Complex Traits Genetics Virtual Lab (CTG-VL) platform (Methods), which revealed significant correlations (p FDR < .05) between IBS types and 278 other traits across the GI, psychiatric, musculoskeletal, respiratory, cardiovascular, and metabolic domains (full results in Supplementary Table 2).Strongest correlations were observed for selfreported IBS and stomach or abdominal pain (Figure 3).Indeed, among different IBS types, IBS-C is the one that most often departed from others in its genetic similarities with other traits, illustrated for instance by reduced or augmented genetic correlations with stomach or abdominal pain, laxatives, and back pain (Figure 3).Among non-GI and non-pain-related traits, the conditions most closely correlating with IBS were psychiatric conditions, including generalized anxiety disorder (r g ¼0.51-0.61;P ¼ 2.2Â10 -2 ), major depressive disorder (r g ¼0.42-0.49;p 1.9Â10 -5 ), and neuroticism score (r g ¼0.35-0.56;P 4.1Â10 -4 ), all 3 showing stronger mean correlations with IBS-M and IBS-D than IBS-C.As an interesting and new observation, significant genetic correlations with IBS (particularly the IBS-M subtype) were noted for cardiovascular diseases (CVD), including family history of heart disease (measured as illness of mother and of siblings; r g ¼0.29-0.45;P 1.3Â10 -3 ) and high blood pressure (BP, measured as illness of siblings; r g ¼0.32-0.38;P 2.2Â10 -2 ), among others (Figure 3, Supplementary Table 2).

Annotation of Genomic Risk Loci and Fine Mapping
Sixty-four genome-wide significant SNP associations (P 5Â10 -8 ) were detected in Rome III IBS GWAS metaanalyses, corresponding to 4 independent loci: 2 in IBS (tagged by rs2035380 and rs9517497 on chromosome 11 and 13), 1 in IBS-M (tagged by rs2048419 on chromosome 8), and 1 common to IBS and IBS-M (tagged by the lead SNP rs6899057 on chromosome 5), as shown in Table 4 and  No genome-wide significant associations were detected for IBS-D or IBS-C (Figure 4).Association signals from the 4 Rome IBS risk loci (regional plots in Figure 5) showed no heterogeneity (Cochran P-Het>0.05) and consistent genetic risk effects (OR direction) across GWAS performed individually in UKBB and LL (Table 4).Notably, rs2035380 (candidate gene BMAL1) and rs2048419 (CLDN23/MFHAS1) represent novel associations identified in this study, respectively, for IBS and the IBS-M subtype.In addition, although genome-wide associations for rs6899057 (ADCY2) and rs9517497 (DOCK9) had been previously reported for IBS, 10 the former is also linked here to IBS-M for the first time.
To gain initial biologic insight, we adopted a phenomewide association study approach to screen Rome III IBS risk loci for known associations (P 5Â10 -8 ) reported in previous GWAS studies (Methods).Similar to the results obtained from genetic correlations, these analyses revealed the 2 novel risk loci, tagged by rs2048419 (CLDN23/ MFHAS1) and rs2035380 (BMAL1), had both been previously associated with multiple traits including cardiovascular (hypertension, coronary artery disease), psychiatric (anxiety-like), and metabolic (body mass index) traits (Figure 6).
A Bayesian approach to fine-mapping was used to try and identify the most likely causative variants from each locus (Methods).This returned 4 independent credible sets of SNPs (range, 20-99 markers each) with high posterior probabilities to include causative variants (posterior probability, 82.0%-91.4%).However, no individual variant was estimated to be causative with high confidence (all variants were lower than 8.0% probability), and the identity of the individual causal variant from each locus remained elusive.Although there was no coding variant among SNPs from the credible set of each locus, several mapped near regulatory regions and transcription factors binding sites (Table 5).

P>.05 P<.05
Figure 2. Analysis of Rome III IBS comorbidities in UKBB.Selected ICD10 diagnoses at the 3-digit code (trait-specific) level are reported, with corresponding OR and 95% confidence interval versus asymptomatic control subjects, derived from logistic regression adjusted for sex and age.Nonsignificant OR (p FDR !.05)are shown as empty circles.Full results are reported in Supplementary Table 1.
Furthermore, many were associated with differential expression (expression quantitative trait loci [eQTL]) of multiple genes in various tissues from the Genotype-tissue Expression (GTEx) Portal (Figure 7; see Methods), including eQTLs in tissues of the GI tract (for rs2048419 and rs9517497 loci) and heart (for rs2048419).

Gene Mapping and Prioritization of Candidate Genes
Positional and cis-eQTL functional mapping (FUMA) and annotation of GWAS (Methods) resulted in a total of 8 protein-coding genes as plausible causative candidates at the 4 IBS and IBS-M risk loci (Table 4).Their expression in different human tissues and cell types are shown in Figures 7 and 8. Previous GWAS associations of these genes with various traits are listed in Supplementary Table 3.
Single genes were mapped to 3 risk loci, hence these represent bona fide candidates to be causative of the observed associations at respective loci: BMAL1 (basic helixloop-helix ARNT Like 1) at the new rs2035380 locus for IBS, ADCY2 (adenylate cyclase 2) at the rs6899057 locus in common for IBS and IBS-M, and DOCK9 (dedicator of cytokinesis 9) at the rs9517497 locus for IBS (Table 1).BMAL1, mapped positionally and via eQTL, codes for a key clock protein that controls circadian rhythms and is highly Figure 3. Genetic correlations (r g ) across Rome III IBS subtypes and other traits.IBS, its subtypes, and other traits are reported based on respective significant (p FDR < .05)genetic correlations, divided into different domains.Empty cells correspond to nonsignificant results or r g between the same traits.Full results are reported in Supplementary Table 2.
expressed in the brain and the GI tract (colon, esophagus, small intestine, and stomach) (Figure 7), and in multiple cell types (Figure 8).This gene has been implicated via GWAS in several other traits, including CVDs and psychiatric disorders (Supplementary Table 3).DrugBank search revealed BMAL1 regulators are targeted by antibacterial drugs with indications for IBS-D (rifaximin, DrugBank ID: DB01220) and traveler's diarrhea (norfloxacin and ofloxacin, Drug-Bank ID: DB01059 and DB01165).For ADCY2, mapped to the rs6899057 locus via eQTLs, considerable expression is observed in brain tissues and cells (Figures 7 and 8) and in myocytes (Figure 8).Previous GWAS studies linked this gene to IBS and other traits, including GI diseases and CVDs, possibly via cyclic adenosine monophosphate signaling pathway (Supplementary Table 3).The encoded enzyme (AC2) is targeted by the drug colforsin indicated for smooth muscle relaxation and vasodilation (DrugBank ID: DB02587).Finally, at the rs9517497 locus, DOCK9 is involved in GTPase activity regulation and brain development, showing high expression in brain and GI tissues (colon, esophagus, small intestine, and stomach) (Figure 7), and muscle and endothelial cell types (Figure 8).It is associated with eQTLs in the esophagus mucosa (Figure 7) and has been previously linked to IBS and many other traits via GWAS (Supplementary Table 3).
Five candidate genes were mapped to the rs2048419 locus for IBS-M, the second novel locus identified in this study (Table 4).We aimed to prioritize the best candidates at this locus based on functional evidence at multiple levels, including eQTL data for GI and brain tissues, gene expression, previous GWAS findings, mechanistic links to IBS, relevant phenotypes from knockout mice, and DrugBank data (Methods).Multiple functional annotations used for gene prioritization in the rs2048419 locus are reported in Table 6.By these means, CLDN23 (claudin 23), which encodes a tight-junction protein that plays a role in gut barrier function, and MFHAS1 (multifunctional ROCO family signaling regulator 1), which modulates immune responses, were deemed the most likely causative genes.Both show expression in the GI tract (Figure 7; including expression in specific colonic cell types, such as enteroendocrine cells, Figure 8), are associated with eQTLs in GI tissues via finemapped credible variants (Table 6), and are linked via GWAS data to other traits, such as neuroticism and hypertension (Supplementary Table 3).

Genetic Link to Heart-Related Traits: Sensitivity and Conditional Analyses
Given the observed associations between IBS and heartrelated traits, we sought to gain further insight from 2 additional, alternative analyses of their relationship.First, using UKBB data, we repeated our IBS Rome III GWAS analyses excluding all participants with diseases of the circulatory system (ICD10 chapter IX; diagnoses I00-I99), leaving 5993 IBS cases and 15,041 control subjects (Methods).This analysis revealed the genetic architecture of Rome III IBS to be virtually identical to that estimated devoid of heartrelated conditions (r g >0.98 and P .4Â10 -14 , across all Rome III IBS subtypes based on LDSC analysis of GWAS data).At the same time, the genetic risk effect estimates (ORs) associated with the 4 GWAS signals remained largely unchanged (despite the reduction in sample size and, expectedly, the strength of the associations) (Table 7).Next, we performed a proper conditional analysis to assess again, based on a different approach, whether the genetic underpinnings of Rome III IBS are independent of heart- disease predisposing factors.For this purpose, using mtCOJO 20 (Methods), we conditioned our Rome III IBS GWAS meta-analysis results on the use of CVD medications (as proxies for their indications), using GWAS summarylevel data from Wu et al 21 (Methods).As shown in Table 8, Rome III IBS GWAS signals were not attenuated on   conditioning, and the genetic architecture of IBS subtypes remained substantially unaltered, as assessed via LDSC analyses of conditioned and unconditioned GWAS data (r g >0.92 and P 7.4Â10 -263 , across all Rome III IBS subtypes).Similarly, conditioning on antidepressants as proxies for anxiety and mood disorders (conditions known to share underlying genetics with IBS) did not result in substantial changes of the genetic risk effects (Table 8; and r g >0.94 and P 1.0Â10 -300 across all Rome III IBS subtypes).Overall, these results suggest that the genetic predisposition to Rome III IBS is not dependent on co-occurring heart-related diseases (and/or anxiety and mood disorders).Instead, the detected GWAS risk signals are intrinsic to the genetic architecture of IBS.

Rome III IBS Polygenic Risk Scores
Finally, we computed and tested Rome III IBS polygenic risk score (PRS), to assess the relevance and translational potential of our GWAS results in predicting IBS risk.For this purpose, in the absence of other large cohorts with Rome III and genotype data available, we used phenotype and genotype data from the independent set of UKBB participants who did not fill the digestive health questionnaire (DHQ nonresponders), and who were therefore not included in the Rome III IBS GWAS (Figure 1).In this dataset, IBS cases were identified based on hospital-inpatient records and self-reported diagnoses, whereas control subjects were the remainder of DHQ nonresponders (Methods).Two similarly sized nonoverlapping fractions (n ¼ 146,771 and 146,772) of the DHQ nonresponder dataset were then used, respectively, for computation and testing of Rome III IBS PRS (see Table 9 for demographics).
Using the PRSice-2 22 pipeline (Methods), we computed a PRS model based on the inclusion of 51,509 SNP markers that were best able to differentiate IBS cases (n ¼ 5335) from control subjects (n ¼ 141,436) in the first dataset (P ¼ 6.0Â10 -13 ).Using this marker set in the testing group (5336 IBS cases and 141,436 control subjects), we detected significantly higher PRS values in IBS cases compared with control subjects (P ¼ 7.9Â10 -22 ).Notably, IBS prevalence progressively increased across the PRS distribution, reaching 4.5% in the top 1% percentile (OR, 1.34; P ¼ 1.1Â10 -3 vs the rest of the population) (Figure 9).

Discussion
We report a large-scale, population-based survey of clinical and genetic data in relation to IBS and its subtypes according to the gold standard Rome III criteria.Based on GWAS meta-analyses of 2 independent cohorts comprising 101,884 individuals of European ancestry, we show that Rome III IBS and its subtypes are heritable traits, with SNPbased heritability ranging between 10% and 13%.This is higher than 5.8% SNP heritability reported in previous larger IBS GWAS meta-analyses of multiple heterogeneous IBS definitions, 15 suggesting that standardized Rome III criteria might better capture the contribution of core common genetic risk factors to IBS.The incomplete overlap of GWAS associations signals detected here and in previous studies may also indicate that different definitions of IBS (ICD10 diagnoses from electronic medical records, selfreported doctor diagnoses, Rome criteria, or all combined) may capture slightly different genetic underpinnings. 9At the same time, the genetic architecture of Rome III IBS is largely shared across subtypes, including those at the opposite end of the stool consistency spectrum (as shown in genetic correlations of IBS-D and IBS-C).Rome III IBS also shows considerable genetic overlap with other GI, psychiatric, and extraintestinal pain-related conditions that, as shown in previous studies, 23,24 are often comorbid as in our UKBB ICD10 analyses.Novel from this study is the observation of strong genetic correlation between IBS and heart-related conditions and traits, including family history of heart disease and hypertension, ischemic (coronary) heart disease, and angina pectoris, ranging between 20% and 45%.This novel link between IBS and cardiovascular traits was further consolidated using 2 independent, different analytical approaches: GWAS devoid of individuals with heart diseases, and GWAS conditioned on CVD medications (as proxy for CVD) both provided evidence that the genetic risk effects observed for Rome III IBS are not dependent on heart-related comorbid conditions, nor on their genetic underpinnings.Showing correlation with all Rome III subtypes, the genetic (sibling history of) hypertension was most pronouncedly similar to that of IBS-C.[27] This was evident across a series of heart-related conditions including hypertension, venous thromboembolism, angina, peripheral artery disease, myocardial infarction, stroke, and arrhythmias.In relation to the latter, it is worth noting that genetic commonalities with IBS have also been described for the SCN5A gene (coding for the NaV1.5 sodium channel with a pacemaker role in myocytes and gut interstitial cells of Cajal), whose rare defective variants play a role both in Brugada syndrome and IBS-C. 11,28In addition, other crossmechanisms may link IBS to heart-related traits, including alterations of the autonomic/sympathetic control of BP and cardiovascular activity, as observed in patients with IBS more severe symptoms. 29Therefore, similar to what has been previously shown for anxiety and mood disorders, pleiotropic mechanisms seem to be important also for the predisposition to IBS and CVDs, something that is highlighted here for the first time.
Evidence linking Rome III IBS to cardiovascular traits also came from the analysis of individual GWAS signals.1][32][33] Although multiple genes were mapped to the rs2048419 locus, only 1, BMAL1, is associated with rs2035380.BMAL1 (or ARNTL) codes for a core transcriptional factor that controls circadian rhythm and chronophysiologic body mechanisms, including digestive and cardiac functions. 34,35Of pharmacologic interest, BMAL1 regulators are targeted by antibiotic drugs used to treat IBS-D and traveler's diarrhea (eg, rifaximin, norfloxacin, and ofloxacin). 6,36BMAL1-deficient mice show altered diurnal oscillations in gut microbiota composition associated with changes in the colonic concentration of short-chain fatty acids and bile acids, microbial metabolites known to be NOTE.Conditional analyses used GWAS summary-level data on the use of CVD and antidepressant medications, obtained from Wu et al 21 (PMID:31015401).CoA, coenzyme A; CVD, cardiovascular diseases; GWAS, genome-wide association study; HMG, 3-hydroxy-3-methylglutaryl; IBS, irritable bowel syndrome; IBS-M, irritable bowel syndrome mixed subtype; OR, odds ratio (with 95% confidence interval) of the associated lead SNP for the discovery IBS and IBS-M GWAS; OR-cond, odds ratio (with 95% confidence interval) of the associated lead SNP after conditioning IBS/IBS-M GWAS results on medication use.
8][39] In line with this, functional studies demonstrated GI clock disruption is associated with abnormal colonic motility. 40Altogether, these findings suggest that similar pathogenetic mechanisms may contribute to Rome III IBS and heart disease via genetic disruption of the circadian system, warranting further investigations.
Interesting observations also came from the functional annotation of the other 3 GWAS risk loci.Association signals at rs6899057 on chromosome 5 and rs9517497 on chromosome 13 had been previously reported, 15 and correspond to single gene annotations (respectively ADCY2 and DOCK9).In the context of findings from this study, it is noteworthy that ADCY2 (adenylyl cyclase isoform II) is the target of colforsin (NKH477), a compound that has smooth muscle relaxant and vasodilation properties with indication for acute heart failure. 41,42Moreover, ADCY2 mutations have been associated with congenital heart disease and abnormal calcium signaling. 43Of relevance to IBS, this gene acts downstream of G-protein coupled receptors, such as the muscarinic acetylcholine 44 and serotonin (5-hydroxytryptamine) 45 receptors, which play key roles in the autonomic control of GI motility 46,47 and cardiovascular function, 48,49 which are pathways targeted by current drugs used to treat IBS. 6Hence, AC2 may hold potential for IBS therapeutic exploitation.
At the novel rs2048419 IBS-M locus on chromosome 8 CLDN23 and MFHAS1 were prioritized among 5 candidate genes mapped to the region.CLDN23 encodes the nonclassical tight-junction claudin-23 protein that has been recently demonstrated, both in vitro and in vivo, to modulate the intestinal epithelial barrier. 50This occurs via interaction of claudin-23 with other claudin isoforms (particularly claudin-3 and -4) at the luminal surface of intestinal epithelial cells, leading to regulation of paracellular permeability to ions and macromolecules and strengthening of barrier function. 50Changes in the expression of other claudins have been documented in patients with IBS (especially IBS-D and IBS-M), and linked to severity of GI symptoms. 51,52Hence, it is possible that similar (still uncharacterized) genotype-drive mechanisms are responsible for the putative involvement of claudin-23 in IBS.4][55] The latter drives inflammation on recognition of microbial components (eg, lipopolysaccharide), and multiple studies report increased toll-like receptor expression in the intestinal mucosa of patients with IBS. 56,57The interplay between gut barrier disruption and immune overactivation is hypothesized to be linked to IBS progression via sensory afferent overstimulation, resulting in pain dysfunctions. 51Interestingly, CLDN23 and MFHAS1 are expressed in intestinal enteroendocrine cells, which are key sensors of the intestinal environment and play important roles in modulating gut-brain sensory and immune responses. 58Additional genetic and functional studies are  needed to unequivocally establish weather CLDN23 or MFHAS1 or both contribute to IBS pathophysiology.Rome III IBS GWAS data were also tested for their potential to contribute to IBS risk prediction, via PRS computation and testing.Despite the lack of accurate Rome III classifications in the independent UKBB dataset used for this purpose, our PRS model still allowed the identification of individuals at higher risk of IBS (self-reported or from hospital admissions), with those in the top 1% of PRS distribution being exposed to more than 1.3 times higher risk of disease.Similar results were recently obtained using stool frequency PRS to capture IBS in UKBB, 14 which collectively indicates that, once refined and integrated with other predictive markers, IBS PRS may find translational scope for clinical utility: they may contribute to algorithms designed for differential diagnosing GI diseases (eg, IBS vs inflammatory bowel disease) and/or the decision-making process associated with costly and/or invasive procedures (endoscopy), with positive repercussions on the health care system and the patient.
Finally, although subtype-specific analyses did not reveal major differences, we observed IBS-C was the trait that most often departed from other subtypes in its genetic correlations with other conditions, suggesting specific genetic risk factors and mechanisms may be associated with constipation.Nevertheless, no GWAS significant signals were detected for Rome III IBS-C and IBS-D subtypes: it cannot be excluded that this is caused, for instance, by the current sample size that is insufficient (ie, smaller than IBS-M) to capture weak individual risk effects.Indeed, the measured heritability of these traits seems to be slightly lower than that of IBS-M.This issue can be better addressed in follow-up, possibly larger GWAS meta-analyses in independent cohorts, or via alternative gene-hunting strategies based on the analysis of endophenotypes, such as stool frequency and consistency. 14mportant study limitations need to be acknowledged, including: (1) having studied IBS defined according to the older version of the Rome III Criteria, which warrants evaluation of current findings in patients defined according to stricter Rome IV (or possibly upcoming Rome V) criteria when suitable large cohorts with these data become available; (2) GWAS analyses not performed for ancestries other than European; (3) minimal understanding of the pleiotropic mechanisms linking IBS to cardiovascular traits; (4) sample size still likely inadequate to fully capture genetic liability for IBS subtypes; and (5) PRS testing performed on suboptimal IBS classifications because of the lack of other suitable Rome III IBS cohorts with genetic data available.
In summary, we show that IBS types defined according to Rome III criteria are heritable traits whose genetic architecture is shared with other polygenic traits, including CVDs highlighted here for the first time.We also identify novel candidate genes and plausible mechanisms that represent potential actionable targets warranting further investigation in follow-up studies.

Materials and Methods
The study design is summarized in Figure 1.Detailed methodologic description is provided below.
Study Cohorts and Definition of Rome III IBS Cases and Control Subjects UK Biobank.The UKBB is a population-based longitudinal cohort of >500,000 UK individuals (aged 40-69 years) with available genotype, demographic, and health-related data. 59n exclusion of individuals of non-White ethnic background (data field 21000), sex mismatches (data fields 31 and 22001), and low-quality genotypes (data field 22010), we focused on data from the DHQ (which includes a Rome III module for IBS), available for a subset of UKBB participants who responded to the invitation to participate in this followup.After exclusion of all individuals with other potentially confounding GI diagnoses that could result in IBS-like symptoms, including celiac disease, inflammatory bowel disease (and others listed in Lifelines.LL is a multidisciplinary prospective populationbased cohort study examining in a unique 3-generation design the health and health-related behaviors of 167,729 persons living in the North of the Netherlands. 60It uses a broad range of investigative procedures in assessing the biomedical, sociodemographic, behavioral, physical, and psychological factors that contribute to the health and disease of the general population, with a special focus on multimorbidity and complex genetics.In this study, we used health-related and genotype data for IBS GWAS analyses in an LL subset with ethical approval for genetic studies.Similar to UKBB, 1990 IBS cases (797 IBS-M, 561 IBS-D, and 475 IBS-C) and 10,518 asymptomatic control subjects were identified according to the Rome III criteria from digestive questionnaire data.The LL protocol was approved by the UMCG Medical ethical committee with protocol number 2007/152.

Comorbidity Analysis in UKBB
For comorbidity analysis, we analyzed 250 UKBB ICD10related data-fields (from 41202-0.0/0.65 to 41204-0.0/0.183) across 89,376 individuals, totaling 2,234,400 datapoints.We tested 1116 ICD10 traits, coded for diseases (chapters I to XIV) and laboratory findings (chapter XVIII), for their differential risk across Rome III IBS cases and asymptomatic control subjects.Associations for the number of recorded ICD10 comorbidities (excluding K58 for IBS) were tested in a linear regression model, adjusting for sex and age in R v3.6.2 (https://www.r-project.org/).ORs were calculated based on logistic regression models (also adjusted for age and sex) at 2 ICD10 levels: chapters (hierarchically subdivided into domains; see https://icd.who.int/browse10/2019/en#); and 3-digit code level (including only 96 conditions with !1% of prevalence in the DHQ subset of UKBB).To control for type I error, the Benjamini-Hochberg false discovery rate (FDR) method was applied in the calculation of statistical significance for multiple tests correction (corrected a¼0.05).

Genotyping, Quality Control, and Genotype Imputation
Similar quality control pipelines and imputation methods were applied to both UKBB and LL cohorts, and are described as follows.
UK Biobank.Genome-wide genotyping of UKBB participants was carried out using custom Axiom (Affymetrix UK Biobank and UK BiLEVE) arrays.Nongenotyped variants were imputed centrally by UKBB researchers, using Haplotype Reference Consortium, UK10K, and 1000 Genomes (1KG) as reference panels.Genetic quality control procedures and PC analysis performed by UKBB researchers are thoroughly described by Bycroft et al. 59 Quality control included checking for deviations from Hardy-Weinberg equilibrium, batch and plate effects, sex effects, and array effects across control replicates.Only SNPs with a call rate >0.99, imputation quality score INFO !0.9, minor allele frequency !0.01, nonindels, nonmultiallelic, and in autosomal chromosomes were tested in association analyses, adjusted for the top 10 PCs (data field 22009).Only Lifelines.Genome-wide genotyping of LL participants was carried out within the UMCG Genotyping Lifelines Initiative.A total of 38,020 participants were genotyped using Infinium Global Screening Array Multi-ethnic Disease Version 1.0 (Illumina). 60Quality controls were based on checking Hardy-Weinberg equilibrium and call rate >0.95 using PLINK 1.9, sex and expected relationships from pedigree, and PC analysis to check for population outliers (as described in http://wiki-lifelines.web.rug.nl/doku.php?id¼ugli).Before imputation, genotypes were prephased using SHAPEIT2 61 and aligned to the reference panels using Genotype Harmonizer 62 to resolve strand issues.Imputation was performed with Beagle 3.1.0 63and Minimac, 64 based on whole-genome sequencing data from the Haplotype Reference Consortium reference panel.Only SNPs with INFO !0.9, minor allele frequency !0.01, nonindels, nonmultiallelic, and in autosomal chromosomes were tested in 1990 samples identified as IBS cases and 10,518 asymptomatic control subjects.

Rome III IBS GWAS and Meta-Analyses
Rome III IBS (and subtypes) GWAS analyses were performed with a logistic mixed model using SAIGE v0.42.1 65 adjusted for age, sex, 10 top PCs, and array (only for UKBB analyses).Individual IBS subtype GWAS summary statistics were inspected with the R package EasyQC v9.0, 66 to check for data integrity, remove invalid or unmapped markers, and harmonize SNP identifiers (rs IDs) and allele strand coding across datasets.Markers with allele mismatches or with allele frequency deviating >0.2 compared with the Haplotype Reference Consortium reference panel were excluded.IBS-subtype specific meta-analyses were performed using a fixedeffect inverse-variance-based analysis with METAL v2011-03-25, 67 including 6,311,313 markers in common between UKBB and LL individual GWAS.A Cochran Q test was performed to assess heterogeneity of effect-size estimates between individual GWAS, and markers with a P-Het <0.05 were excluded from subsequent subtypespecific analyses.

Bayesian Fine-Mapping Analysis
A Bayesian fine-mapping analysis was carried out using FINEMAP v1.4 68 , to determine a minimum set of variants (credible set) containing likely causal variants for each associated risk locus.Posterior inclusion probability was estimated for each credible variant.As input, all genetic variants located within each risk locus were extracted, and the local linkage disequilibrium (LD) structure was calculated using genotypes from European UKBB individuals as reference.

Annotation of GWAS Risk Loci and Gene Mapping
For annotation of GWAS risk loci and gene mapping, we used the online tool FUMA v1.3.5 69 (https://fuma.ctglab.nl/tutorial).All loci at genome-wide significance (P 5Â10 -8 ) were annotated by FUMA, with lead SNPs being defined as those with the lowest P values in each genomic risk locus.Manhattan plots were produced using the R package qqman 70 and Locuszoom 71 was used for regional visualization of risk loci.Independent significant SNPs were identified based on LD with the lead SNPs (R 2 >0.6, using 1KG-EUR as reference panel) and P value (P 5Â10 -8 ).The maximum distance between LD blocks to merge into a locus was set to 250 kb.Gene mapping at each risk locus was based on position (distance <10 kb) and cis eQTL (p FDR < .05)analyses, excluding the major histocompatibility complex region, in view of its LD complexity.The analyzed eQTL database included 54 tissues from GTEx v8 72 (https://gtexportal.org/).

SNP Heritability and Genetic Correlation Analyses Via LDSC
SNP-based heritability (h 2 SNP ) and genetic correlations (r g ) with other 1400 traits were estimated using LDSC, 73 as implemented in the CTG-VL 74 platform.Traits with <1000 cases, <5000 overall participants, and duplicated were removed.P values were adjusted for multiple comparisons using the FDR method (corrected a¼0.05).

Phenome-Wide Association Study
For each IBS risk locus, lead SNPs and/or their LD proxies (R 2 >0.8) were used to screen for previously reported GWAS signals (P 5Â10 -8 ), inspecting data from OpenTargets, 75 PhenoScanner v2, 32 and GWAS ATLAS. 33omains associated with individual GWAS findings were harmonized according to GWAS ATLAS classifications.

Prioritization of Candidate Genes
For the novel risk locus with multiple candidate genes (rs2048419), we attempted to identify those most likely to play a causative role based on functional evidence from finemapping (FINEMAP 68 ), eQTL and gene expression (GTEx 72 and Human Protein Atlas 76 ), GWAS databases (Open-Targets, 75 PhenoScanner v2, 32 and GWAS ATLAS 33 ), Drug-Bank (https://go.drugbank.com/),and literature mining.Genes were prioritized based on the following criteria: (1) considerable expression (set as log 2 transcript per million !3) in GTEx tissues of interest (GI and/or brain), (2) credible SNP-eQTL colocalization in tissues of disorder of gut-brain interaction interest (GI and/or brain), (3) known association with relevant traits from previous GWAS studies, (4) plausible mechanistic link to IBS pathophysiology; and (5) relevant phenotype from knockout/knockdown models or human studies.

Rome III IBS Risk Loci Sensitivity Analysis of CVD Traits
Reevaluation of Rome III IBS loci association signals in UKBB devoid of heart-related conditions was obtained by removing from the analyses individuals with diagnoses from the ICD10 chapter IX diseases of the circulatory system.These included: I10 essential (primary) hypertension, I20 angina pectoris, I25 chronic ischemic heart disease, I44 atrioventricular and left bundle-branch block, I48 atrial fibrillation and flutter, I73 other peripheral vascular diseases, I83 varicose veins of lower extremities, I84 hemorrhoids, and I95 hypotension.As previously, association testing was carried out with SAIGE v0.42.1, 65 adjusted for age, sex, 10 top PCs, and genotyping array.In total, 15,041 individuals with chapter IX diagnoses were removed from asymptomatic control subjects, whereas 5993 were removed from Rome III IBS cases (5993 from IBS, 3380 from IBS-M, 1475 from IBS-D, and 998 from IBS-C).

Multitrait Conditional Analysis
Multitrait-based Conditional and Joint analysis (mtCOJO 1.93.2 beta 20 ) was used to generate Rome III IBS GWAS summary statistics (outcome) conditioned on GWAS data for the use of CVD and antidepressant medications (exposures).mtCOJO estimates the effect of the exposure on the outcome either by generalized summary-data-based mendelian randomization or from genetic correlation (rg) analysis, when there are not enough LD independent (r 2 <0.05) genome-wide significant SNPs for the exposure or outcome traits (at least 10 required by default).mtCOJO analyses were carried out on IBS and IBS-M GWAS meta-analyses data versus GWAS of CVD medication traits derived from Wu et al, 21 including: C01D vasodilators used in cardiac diseases, C02 antihypertensives, C03 diuretics, C07 beta blocking agents, C08 calcium channel blockers, C09 agents acting on the renin-angiotensin system, C10AA 3-hydroxy-3methylglutaryl coenzyme A (HMG CoA) reductase inhibitors, and N06A antidepressants, according to the drug coding classification of the World Health Organization Anatomical Therapeutic Chemical Classification System.

Polygenic Risk Score Analysis
We adopted a 2-step approach to PRS computation and testing in independent subsets of UKBB, only studying participants who did not fill the DHQ questionnaire (DHQ nonresponders) and were, thus, not included in the original Rome III GWAS meta-analyses.In the absence of Rome III data for these individuals (Rome III modules are part of the DHQ questionnaire), IBS cases were identified based on hospital-inpatient records (ICD10 diagnosis K58; data fields 41202 and 41204) and self-reported doctors' IBS diagnosis (1154 code; data field 20002), whereas control subjects were the remainder of QCed DHQ nonresponders.Cases and control subjects were randomly assigned to 2 independent and similarly sized groups for PRS computation (5335 cases and 141,436 control subjects) and testing (5336 cases and 141,436 control subjects).
For PRS computation, according to the standard LD pruning and P value thresholding of PRSice-2 v2.2.11.50 22 pipeline, varying numbers of SNPs with GWAS P values ranging from 5Â10 -8 to 1 were tested using the Rome III IBS GWAS summary statistics as the base file.Using default PRSice-2 pruning parameters (clump-kb¼250, clump-p¼1, clump-r 2 ¼ 0.1) and 1KG-EUR LD as the reference panel, a total of 214,641 LD pruned variants entered the analyses, after exclusion of poor-quality markers (genotype missingness geno <0.1 and imputation quality INFO<0.9) and rare markers (minor allele frequency <0.01).The best Rome IBS PRS model was derived in the computation dataset, by selecting the largest and most significant Nagelkerke R 2 value obtained in the P value thresholding strategy (P(R 2 ) ¼ 6.0Â10 -13 ; GWAS P threshold ¼ .1056),which included 51,509 SNPs.
For the testing group, PRSs were then calculated for IBS and control subjects as the sum of effect sizes of the 51,509 markers that fitted the best Rome IBS PRS model.PRS distributions were scaled (mean, 0; standard deviation, 1) using the "scale"' R function.A multivariate linear regression (adjusted for sex, the first 10 top genetic PCs, and UKBB genotyping array) was used to determine the significance of the difference between the mean Rome IBS PRS in IBS cases and control subjects.Individuals within a given magnitude of increased PRS in top percentiles (binned using the "ntile" dplyr R function) were compared in relation to the rest of the population for their risk to IBS in a logistic regression including previous covariates, and FDR was applied in the calculation of statistical significance for multiple tests correction.

Figure 1 .
Figure1.Infographic summarizing the derivation of study cohorts (IBS cases and control subjects) and their inclusion in different analyses.Digestive health questionnaire (DHQ) responders are described as "DHQþ" and nonresponders as "DHQ-."

Figure 4 .
Figure 4.No genome-wide significant associations were detected for IBS-D or IBS-C (Figure4).Association signals from the 4 Rome IBS risk loci (regional plots in Figure5) showed no heterogeneity (Cochran P-Het>0.05) and consistent genetic risk effects (OR direction) across GWAS performed individually in UKBB and LL (Table4).Notably, rs2035380 (candidate gene BMAL1) and rs2048419 (CLDN23/MFHAS1) represent novel associations identified in this study, respectively, for IBS and the IBS-M subtype.In addition, although genome-wide associations for rs6899057 (ADCY2) and rs9517497 (DOCK9) had been previously reported for IBS,10 the former is also linked here to IBS-M for the first time.To gain initial biologic insight, we adopted a phenomewide association study approach to screen Rome III IBS risk loci for known associations (P 5Â10 -8 ) reported in previous GWAS studies (Methods).Similar to the results obtained from genetic correlations, these analyses revealed

Figure 4 .
Figure 4. Manhattan plots of subtype-specific IBS UGWAS meta-analyses.Each circle denotes a marker with regard to its physical location (based on the Genome Reference Consortium Human Build 37) and associated -log 10 P value.SNPs reaching the genome-wide significance threshold (P ¼ 5Â10 -8 , indicated with a horizontal red line) are colored in green.Suggestive significance threshold (P ¼ 5Â10 -6 ) is indicated with a horizontal blue line.The lead SNP rs ID and its respective best candidate protein-coding gene symbol are annotated in green.

Figure 5 .
Figure 5. Regional plots of associated IBS and IBS-M risk loci.SNPs (dots) are reported with their physical location (based on the Genome Reference Consortium Human Build 37), associated -log 10 P value, and degree of linkage disequilibrium (r 2 ) with the lead SNP (purple diamond labelled with the rs ID).The red line indicates the genome-wide threshold at P ¼ 5Â10 -8 and the blue line indicates the suggestive threshold at P ¼ 5Â10 -6 .Purple lines labelled as "EUR" indicate the genomic location of fine-mapped credible variants.

Figure 9 .
Figure 9. Rome III IBS PRS testing analysis.Prevalence of IBS across PRS percentiles (polynomial regression fitting line in black), and OR relative to various upper percentiles compared with the rest of the population (estimated from logistic regression adjusted for sex, 10 top principal components, and genotyping array).The testing cohort consisted of IBS cases (n ¼ 5336) and control subjects (n ¼ 141,436) identified based on ICD10 K58 (for IBS) and self-reported IBS diagnosis.

Table 2 .
Number of ICD10 Comorbidities in UKBB Subjects

Table 4 .
Rome III IBS Risk Loci Identified Via GWAS Meta-analyses x BP, base-pair position (genome build hg19); CHR, chromosome; CI, confidence interval; EA, effect allele; EAF, effect allele frequency; GWAS, genome-wide association study; IBS, irritable bowel syndrome; IBS-M, irritable bowel syndrome mixed subtype; LL, Lifelines; OA, other allele; OR, odds ratio; SNP, single-nucleotide polymorphism; UKBB, UK Biobank.OR in individual cohorts: magnitude of mean effects described to the risk allele in UKBB and LL individual GWAS; P-Het: Cochran's Q statistics P value for heterogeneity across individual GWAS; Number of credible variants: number of markers in each fine-mapped credible set per risk locus; Gene mapping: genes mapped to the locus, based on FUMA positional ( p ) or cis-eQTL ( e ) analyses.The best candidate genes at each locus are highlighted in bold.

Table 5 .
List of Fine-Mapped Credible Variants That Are Near Regulatory Sites Based on Ensembl Variation Annotation IBS 1.06 (1.04-1.08)6.3E-07

Table 5 .
Continued NOTE.Functional consequence of variants was performed using the Ensembl Variant Effect Predictor tool (https://grch37.ensembl.org/info/docs/tools/vep). OR: mean effects described to the risk allele in the meta-analysis.BP, base-pair position (genome build hg19); CHR, chromosome; CI, confidence interval; EA, effect allele; EAF, effect allele frequency; GWAS, genome-wide association study; IBS, irritable bowel syndrome; IBS-M, irritable bowel syndrome mixed subtype; LL, Lifelines; OA, other allele; OR, odds ratio; SNP, single-nucleotide polymorphism; UKBB, UK Biobank.Figure7.Heat map of gene expression from Rome III IBS risk loci, based on GTEx data.IBS candidate genes from IBS risk loci (locus membership annotated at the top) are reported with their level of expression (mRNA) in 54 human tissues using data from GTEx v8 (https://gtexportal.org).The log 2 (TPMþ1), where TPM (transcripts per million) is the averaged expression per tissue and per gene, is indicated on a color-coded scale.eQTLs (p FDR < .05)linked to lead SNP or variants in high LD (r 2 >0.8) are boxed.

Table 6 .
Functional Evidence Used for Prioritization of rs2048419 Candidate Genes NOTE.Functional evidence used for gene prioritization included: (1) GI/brain expression: high expression (defined as log2TPM!3) in GI and/or brain tissues based on GTEx data; (2) GI/brain eQTL: colocalization between a credible variant and an eQTL in a GI and/or brain tissue; (3) Gene-trait association: association with IBS or other traits in previous GWAS studies via credible variants in the gene region; (4) Mechanistic link: plausible mechanistic link between the gene and IBS; (5) Phenotype from KO model or human studies: relevant phenotype from experimental knockout/knockdown models or human studies.Obs: No genes in this locus were directly targeted by current drugs with IBS-relevant indications based on DrugBank search.eQTL, expression quantitative trait loci; GI, gastrointestinal; GTEx, genotype-tissue expression; GWAS, genome-wide association study; IBS, inflammatory bowel syndrome.

Table 7 .
Sensitivity GWAS Analysis After Removing Individuals with CVD Diagnoses in UKBB aEstimates from Rome III IBS GWAS meta-analyses (UKBB and LL).

Table 8 .
Conditional GWAS Analysis on the Use of CVD and Antidepressant Medications

Table 9 .
Demographics of UKBB Study Cohorts Included in PRS Analyses

Table 10 )
, patients with IBS (N ¼ 22,745 total; 12,335 IBS-M; 5920 IBS-D; 3855 IBS-C) were identified using Rome III DHQ questions (data fields 21025, 21027, 21028, 21031-21034).Asymptomatic control subjects (n ¼ 66,631) were defined as the remainder of the population without abdominal pain/discomfort in the previous 3 months (data field 21025) and without IBS, based on IBS diagnoses from hospital inpatients records (ICD10 K58 code; data fields 41202 and 41204) and selfreported IBS from questionnaire data (data fields 21024 and 20002).The study was approved by the Monash University Institute Research Ethics Committee with protocol number 20326 and performed using UKBB data accessed under application number 17435.

Table 10 .
Exclusion Criteria Applied for UKBB IBS Cases and Control Subjects Data-field Trait individuals with White European ethnic background (categorized as "British," "Irish," or "Any other White background"; data field 21000) and without poor heterozygosity and missing rates (data field 22010) were included, leaving 22,745 Rome III IBS cases and 66,631 asymptomatic control subjects for GWAS analyses.